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ABSTRACT 



The Delay Discriminator, an Optimum Tracking Device, is under current 
study as a new space tracking system. The transient analysis of this 
system requires the solution to a non-linear second-order differential 
equation. The Delay-Lock principle is first explained and the implimenta- 
tion of this principle is described. Analysis of the system is discussed 
and a method of digital phase-plane solution is shown. Computer Flow 
Diagrams and Fortran Source Programs are fully described. With out the 
aid of a high speed digital computer, the analysis of this system would 
be difficult, tedious, and inaccurate. 
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1 . Introduction « 



The advent of the space age has created many problems concerned with 
position measurements and tracking. In this situation, it is sometimes 
necessary to measure the time delay envolved in sending and receiving a 
given signal when the target is at great distance and moving at high 
velocity. The Digital Delay-Lock Discriminator described in this paper 
is a statistically optimum system which measures the delay between two 
correlated waveforms and satisfies the above distance and velocity re- 
quirements; it mathematically operates in a similar fashion to that of 

/3/ 

the Phase-Lock System . The Discriminator is a non-linear system in- 
corporating feedback, and attempts to correlate the delay of the return 
signal with an estimate of the delay generated in the receiver. The 
Delay-Lock is mathematically described by a second-order non-linear 
differential equation, whose solution is not readily obtainable by analog 
methods. The transient analysis of the Delay-Lock Discriminator is the 
main objective of this paper and will be described fully in a later 
section. 

The following section will present a description of the basic 
principles of the Delay-Lock System as first described by Spilker^^^; 
and detail a digital implimentat ion of this system currently operational 
at the Lockheed Missle and Space Company’s Research Labs. This is 
followed by an explanation of phase-plane analysis and methods of graphical- 
ly sketching the transient behavior, 

A further section is devoted to a method of dynamic or sequential 
programming, which allows the computer to step off the operating tra- 
jectories in the phase plane. A CDC 1604 High Speed Digital Computer 
using Fortran Programming provided the transient data which was later 
plotted. 
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This paper then describes a new tracking system applicable to 
many space problems, and presents the transient analysis of its tracking 
operation. 



2, System Operation o 

This section is concerned with the basic theory of operation of 



tion of the basic principles of the Delay-Lock System will first be 
presented. After this, a method of generating a pseudo-random sequence 
to represent the signal is described. In conjunction with this signal is 
the techniques of autocorrelation to obtain the desired discriminator 
characteristic. The mathematics of the autocorrelation procedure is 
then completely detailed. All of the above topics are then encompassed 
in a description of the operating Digital Delay-Lock Discriminator. 

2-1 Theory of the DELAY-LOCK DISCRIMINATOR. 

The Block Diagram of the Delay-Lock Discriminator is shown in Fig. 
2-1. The system is a non-linear discriminator incorporating a multi- 
plier, a low pass filter, and a variable delay element. This delay 
element may be ultrasonic, fixed ferrite cores with multiple taps, or 
as in the case that follows, a shift register operating at variable clock 
rates. 
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Fig. 2-1 Block Diagram of the Delay-Lock Discriminator. 
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In an unpublished paper^ , Spilker has shown that this discriminator, 
or a slightly modified version is a optimum system in that it provides 
the maximum likelihood (a posteriori, most probable) estimate of delay. 

Defining the delay error such that 
where estimate of Xft) and expanding the delayed signal by 

Taylor's Expansion produces the following expression: 



SCt+T)= ) + € 5’ it + ^ • ' 

Zj 

With S (t) normalized to have unity power it can be shown that the output 
of the multiplier X(t) can be given by: 

“ff" = f\Td 6 (ir) + ne C'^) 

where ^ error correcting term and noise term. 

In viewing the above equation, the basic principle of the discrimina- 
tor tracking operation is apparent. Consider a sudden increase in the de- 
lay time T(t1 . Such an increase will then increase the initially small 
delay error ^ a) . As 6(0 increases, the correcting term of the 
multiplier output xU) increases, and thus the estimate of delay f(t) will 
increase and tend to track the input delay. 

A partially linearized network in Fig. 2-2 describes the multiplier 

f(p) 

output. The close loop transfer function l4 f ^ ^ . i. ■ 
is linearized when ^ is constant and describes the locked on region 
when id) is small. 
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Fig. 2-2 Partially linearized circuit for delay-lock system. 



In an unpublished paper^^\ Spilker has shown that this discriminator, 
or a slightly modified version is a optimum system in that it provides 
the maximum likelihood (a posteriori, most probable) estimate of delay. 
Defining the delay error such that £lC6j ^ T( t) 

where estimate of XCt) and expanding the delayed signal by 

Taylor's Expansion produces the following expression: 



S(t+T)= sC<-tf) + e 

With S (t) normalized to have unity power it can be shown that the output 
of the multiplier x(t) can be given by: 

^ f\ Id 6 ( t) + He 

where ^ error correcting term and a noise term. 

In viewing the above equation, the basic principle of the discrimina- 
tor tracking operation is apparent. Consider a sudden increase in the de- 
lay time T(tl . Such an increase will then increase the initially small 
delay error ^ Ct) . As fe«) increases, the correcting term of the 
multiplier output xU) increases, and thus the estimate of delay m will 
increase and tend to track the input delay. 

A partially linearized network in Fig, 2-2 describes the multiplier 

Tf p) 

output. The close loop transfer function l4 f ^ ' - ■ 

is linearized when ^ is constant and describes the locked on region 
when small. 
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Fig. 2-2 Partially linearized circuit for delay-lock system. 
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Spilker shows in his paper that if Stt/is a 
sense) ergodic random variable with zero mean and slowly varing delays 
T(i) and f(i) the loop filter, when optimized, forms the average 
of the multiplier output to obtain: 

E [xct)]-- E([/As(t*T)*nW] i<s‘Cii-f)l 
— K/lRjCT-t) 

Where i?:(t 4) is the derivative of the auto-correlation function of 
s(i) . Hence is not linearly dependent on the delay error but 

functionally dependent upon the error through the differentiated auto- 
correlation function, and thereby causes changes in the effective loop 
gain. 



By considering the discriminator characteristic for a gaussian low- 
pass signal spectrum, it is possible to examine the threshold of interest. 
Such a spectrum and corresponding discriminator characteristic is shown in 
Fig. 2-3. 




Fig. 2-3 Signal power- spectral density and discriminator 

function for Gaussian low-pass spectrum. /!/ 



The discriminator characteristic shows only one lock-on region limited by 
I ^ . As exceeds , the slope becomes negative and 
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produces a decrease in delay estimate ; here the system is unlocked 

and unstable. Since there is only one region of positive slope, there 
can be no ambiquities in the lock on region and as such, does not possess 
the limitation in tracking that a sine wave signal contains. 

The lock-on performance is more difficult to describe, but when the 
system contains a simple low-pass RC filter, the mathematics are described 
by a first-order non-linear differential equation. When neglecting noise 
terms this equation is as follows: 

In a later discussion, a simplified form of this e<]uation will be 
used to describe the transient operation of interest. It is well to 
mention here the loop filter characteristics as proposed by Spilker. The 
equivalent filter circuit is shown in Fig. 2-4. The loop filter has 
shown to be optimum for ramp inputs in the presence of white noise, in 
that is minimizes the total squared transient erroi^ plus the mean squared 
error caused by noise. 




Fis* 2-4 Equivalent loop filter and corresponding 

transfer function. /!/ 



The closed loop transfer function for the system is given by; 
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/~2 Generation of the Pseudo- Random Signal. 

This section determines to a large extent the form of the digital 
Delay-Lock System is based on the work of E?spas^^( 

A signal with pseudo-random characteristics and of maximum length can 
be generated with the use of linear feedback shift registers. However, 
not all feedback shift registers will generate maximum length sequences. 
The characteristics of various feedback arrangements are contained in a 
special set of describing polynomials which are of particular 

concern. '' 

In considering the binary linear sequential network as shown in 
Fig. 2-5, its operating cycle can be expressed by a T matrix. 




Ha If -Adder Operation 

1 + 1=0 
0 + 0 = 0 
1 + 0 = 1 
0+1 = 1 

A*B + A*B 



Fig. 2-5 Binary linear sequential network with 
four delay elements and half-adder. 

States S- , So, S^, S, describe the internal condition of the network as 
the four delay element outputs at any time, ^2** ^3** ^4* describe 

the inputs to the delay elements, and will be the outputs of these elements. 



a unit time delay later. Therefore, the S^’ are Mod-2 (half adder) sums 
of the S^ and the following equations apply: 

Si' = S3 + 

S2* = 
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Expressed in matrix notation, these equations may be represented as 
follows : 
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Or, symbolically, as 

S* = T^S 

where T is called the T-matrix representation of the network. 

Fig. 2-6 shows the state diagram of the present example. The net- 
work generates a maximal length sequence (disregarding the all zero state 

4 

which is trival) representing 2 - 1 = 15 possible states. 




Fig. 2-6 State diagram of network shown in Fig. 2-5. 

It is not necessary to examine point by point the different states of 

the network to determine if any given arrangement will generate maximum 
/4/ 

length sequences. This important characteristic can be determined from 
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the network polynomial 0(x), where 0(x) is defined as follows: 



0(x) = jl - xl j 
or in this present example; 

0 



0(x) = 



-X 

1 

0 

0 



-X 



1 
0 

1 -X 



0 



1 



= X + X + 1 



In order for a network to generate a maximal length sequence, its 

characteristic polynomial 0(x) must be irreducible and not be a divisor 
k ^ n 

of X - 1 for any interger k \ p - 1 where n is the number of delay 
elements. 

The output of any delay element over the complete cycle period 
results in the wave shape shown in Fig, 2-7, 
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Fig. 2-7 Output wave-form from third delay element 
of feedback network shown in Fig, 2-5, 

It is such a sequential network that provides the desired signal 

and estimate of delay time in the Digital Delay-Lock Discriminator, 

This network has a characteristic polynomial 0(x) as follows: 

0(x) = + 1 

This network is shown in Fig, 2-8, and the output from any one delay 
element constitutes a continuous signal s(t) ^ which has pseudo- random 

9 

characteristics and a maximum length sequence M = 2 - 1 or 511 pulses. 
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Fig. 2-8 Maximal length sequential network equivalent to 
linear feedback shift-registers utilized in the 
Digital Delay-Lock Discriminator. 

Although these Sequential Networks have the desirable properties 
for generating a maximal length pseudo-random signal; there is one 
characteristic, the auto-correlation function which has yet to 

be described. The autocorrelation function for the network’s sequential 
output signal will be discussed in the following sub-section and shown 
to be ideal, in that it results in a unique and highly desirable Discrim- 
inator Characteristic, 

2-3 Mathematics of the Digital Delay-Lock Discriminator. 

The basic delay- lock discriminator tracking system, covered in 
Section 2-1, can be digitally implemented to track a pseudo-random binary 
signal generated by linear feedback shift register methods previously 
mentioned. The advantage of this new system is that the delay line mention- 
ed in the basic delay- lock system, with its limited tracking range, is re- 
placed by shift registers operating at variable clock rates. 

As mentioned in Section 2-2, a maximum length (2^ - 1) binary sequence 
generated by the output of any state m of an n stage linear feedback shift 
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register is to be the signal, s^(t). The amplitude of s^(t) will 
range from + 1, with a pulse width 'f producing a sequence period of 

seconds. Additive noise effects will in general be ignored in the 
following discussion. 

A simplified block diagram of the digital delay- lock network is 
shown in Fig. 2-9. The binary sequence plus noise is fed to a limiter 
network which converts the input signals to binary form. 




Fig. 2-9 Block diagram of the digital delay- lock 
discriminator tracking system. 

The output of the limiter circuit, r(t) , is then fed to the primary 
logic circuit which consists of two multipliers (AND gates) and a diode 
resister summing net. The output of the logic circuit, x(t is basically 
tei nary (+1,0, -1) and has an average value which allows the delay-lock 
discriminator to track the received signal. The low pass filter follow- 
ing the logic circuit serves to remove as much of the noise and other 
unwanted signals as possible. The output of the filter, T , is the deriva- 
tive estimate of the delay time, proportional to the radial velocity of 
the target reflecting or returning the transmitted signal. This output 
signal is used to control the clock rate of the receiver shift register. 
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Comparison of the clock rates of the transmitter and receiver shift 
registers give a measure of the radial velocity while the Delay is 
measured directly by comparing the instants when both shift registers go 
through a specific state, i.e., theQ,0,0,0,0, — ,0,oQstate. Every 
seconds a delay measurement is available* 

In general, the theory of operation is very similar to that 
previously mentioned in Section 2-1; that is, the received signal and 

the derivative of the transmitted signal, having been delayed an esti- 

A 

mated time T, where multiplied to produce an error correcting voltage . 

However. In this case, the received signal Is represented by s^(t). and 

the delayed version of it represented by s^(t + / ) - s^(t -/) ; both 

signals are then to be multiplied together. The delayed signal is similar 

dS 

to the expression of the time derivative of a continuous signal, — in 

that : cl S . I <,C-t+r)-sCi-r) 

dt -2 T' 

By defining the delay error ^ (t) to be the difference between the 
True and Estimated delays, 6 (t)= Tit) -Tit) , and neglecting 

noise, the output of the multiplier logic circuit, x(t ) > is defined as 
follows: 

■Z e) = Sm r 4 (rt) [s„ (t - r; - s„ « + ^)J 

Being periodic at intervals = 1/ seconds, the received signal 

can be represented by a Fourier series such that: 

cO 

Sfti (^) - Sn Sin tMot) + /)crt 
n-e 

Corresponding to this expression, the difference +7^) ^ ^ 

can be written as follows: 

oO 

Sh(t» * "2 ^ Sih Coj(ni»,-i) 

hr 0 

where 

Asn Cos (nuipV) 6cfi = Acn StriCnWof) 
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and 



By combining the three above equations, x(t, (£ (t)) is as follows; 

x(.i,Ht)) = t y[t, t(t)] 

where D (t)jis defined as the Discriminator Characteristic such that: 

^ ^ mo 

D (t)] = Dtn Cos (nU), t)) 

n»o 

D jn = ”2 (^sn B<q- A-n ds/>) j Dc^) = -“2. (A$n "Bsn Bcn) 

The term can similarly be expressed, and represents the 

intrinsic noise of the system. 

In considering the Discriminator Characteristic, D( ^ ), it is 

seen that this delay correction term is not explicitly dependent upon 
its functional behavior is determined by fixing ^and averaging x(t). 

Thus it is apparent that 0 ( 4 ) is the expected value of the multiplier 
output such that : 

D(t)-- E [* (i), 6)] - E [ S,„ -Y) -Sr,a i- rj] 

This expected value can be expressed in terms of the autocorrelation 

function R ( ^ ) for the binary sequence s (t) so that D( A ) becomes: 
sm m ^ 



Thus the Discriminator Characteristic is defined in terms of the auto- 
correlation functions* As was previously stated, the maximum length 
sequence has a unique and desirable R, ); this function is shown 
in Fig* 2-10^^^. The Discriminator Characteristic is shown in Fig* 2-11: 




Fig* 2-10 Autocorrelation function for maximal 
length binary sequence* 
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Fig, 2-11 Discriminator characteristic for Digital 

Delay-Lock Systems 

It is apparent that D( ^ ) provides a correcting voltage in a 
linear manner when operating in the region \tur , decreases to zero 

is periodic every seconds. 



2-4 The Digital Delay-Lock Discriminator 

This section describes the operating Digital Delay-Lock Discrimina- 
tor System built under the direction of Dr, James J, Spilker, Jr,, at the 
Lockheed Missle and Space Company Research Laboratory, Palo Alto, Calif, 

As shown in Fig, 2-12, the equipment is rack mounted and requires ex- 
ternal analog and digital readout for data analysis. The following major 
components are incorporated in the basic system: 

(a) Rack Mounted Oscillascope 

(b) Transmitter Output Pannel 

(c) Transmitter Shift Register & Logic Circuits 

(d) Receiver Output Pannel 

(e) Receiver Shift Register & Logic Circuits 

(f) Power Control 6e Voltage Controlled Oscillators 

(g) Loop Filter Inputs and Output, and In- lock 
Tracking Light 

(h) Power Supplies 

The basic logic cards (AND/OR) which comprise the major portion of the 
Transmitter and Receiver sections were scrapped from a defunct Polaris 
Missle check out system and as such do not represent the state of the art 
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Fig Z-\Z Digital Delav^Lock Svstpn\ 



in transistor circuitry; however, they provided an inexpensive method 
for realizing the delay- lock system. With the exception of the power 
supplies, the system is completely transistorized; those circuits not 
available originally were built on blank cards. The basic structure of 
the Linear Feedback Shift Registers and associated circuits is shown in 
Fig. 2-14. The primary difference between the two shift registers is 
that the Transmitter shift register operates with a fixed clock rate while 
the Receiver Shift register is controlled by a variable clock rate. The 
clock pulse is generated by a stable crystal oscillator which is fed to a 
doubler or mixer amplifier. This output is fed to a Schmitt Triger cir- 
cuit and then fed to a blocking oscillator. The output is then fed to a 
standard clock-pulse an:plifier and then to the clock inputs of the shift 
register. Although both Voltage Controlled Oscillator networks have a 
manual control for frequency deviation, the receiver VCO is also adjust- 
able by the voltage from the output of the Loop Filter. 

By allowing the difference in clock rates to be large to insure the 
system would not lock-on and track, the error correcting voltage output 
from the filter will conform closely to the system discriminator character- 
istic. When this voltage was fed to an analog plotter, it produced the 
discriminator characteristic shown in Fig. 2-13. This waveform agrees close- 
ly to that predicted for the digital model, the sawtooth discriminator 
characteristic. 




Fig. 2-13 The observed Discriminator Characteristic 
for the Digital Delay-Lock System 
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Fig. Block Diagram of I'^eo Back Shiff f?fe6isrt 








I 



The outputs of both shift registers were fed to a duel-beam oscilla- 
scope to observe the lock-on and track conditions; Fig. 2-15 shows 
both binary sequences in an'’unlocked” and a locked-on and tracking” 
condition. 
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(b) 



Fig. 2-15 Transmitter and Receiver Binary Sequences; 

(a) in an unlocked condition, 

(b) in a lock-on and tracking condition. 

Fig. 2-16 shows an expanded version of Fig. 2-15 





(b) 



Fig. 2-16 Expanded version of Fig. 2-15. 
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The block diagram of the complete digital system is shown in 

Fig. 2-17. The system is operable in the presence of large noise 

signals and was tracking smoothly with an input Signal to Noise ratio 

of approximately - 30 db. Spilker has shown the theoretical limit to 

/!/ 



be-48db. 
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3 . Phase Plane Analysis , 



This section is concerned with the methods of phase plane analysis 
to obtain the transient behavior of the Delay-Lock Discriminator. The 
first topic to be considered is the basic theory of the phase plane and 
its method of application. Following this will be the development of 
the differential equation which describes the lock-on operation of the 
delay-lock system. The final topic to be considered is the computer 
approach to the phase plane solution for two discriminator characteristics, 
the saw tooth and gaussian cases. 

Appendix A contains the flow diagrams and computer programs written 
in Fortran Source Language for the construction of the phase planes of 
interest. 

3-1. Basic Theory and Application of Phase Plane Analysis. 

A large class of Feedback Control Systems and their operation can be 
described mathematically by a differential equation such as: 

X + Pc^,x)z i-Qcx^z) =0 

where P and Q are functions of a signal and its derivative. In phase plane 
analysis methods, the value of ^ is plotted as a function of with the 
individual curves or trajectories for the variation of X and pc an in- 
dication of time. When examining the transient response of a given system, 

# 

the initial conditions of ^ and are plotted and the trajectory (which 
is unique) through this point is the response of the system for these in- 
itial conditions. If the trajectory path converges to a singular point in 
the phase plane, and the system has come to rest with a final value the 

coordinates of the singular point. 

16 / 

Truxel notes three fundamental restrictions on the use of phase plane 
methods : 
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!• "The phase-plane is useful for the analysis of second-order 
systems only " 

2. "The phase-plane can be useful to study only the transient 
performance of a system subject to initial conditions but other wise un- 
excited. " 



3. "The third basic restriction relates to the admissible 
types of non-linearities. The coefficients of x-dot and x can be func- 
tions of x and x-dot, but not of time explicitly." 

Phase plane methods are probably best explained by use of a simple 
example. Consider a second-order differential equation such as: 

X fUJnX = 0 

The time operator t can be eliminated by expressing variations in 
in terms of ^ . Define a new variable^ = ^ and the equation becomes 

Dividing by glvest ^ ^ ^ 2L - - 



now define a slopeO^ such that ^ again substituting gi 



Lves 



Y 



o< 



The phase plane solution to this equation is found by plotting the iso- 
76 / 

dines' or lines of constant slope. This is shown in Fig. 3-1 along 




22 








.w* 



/ 



m.. 



.* 



1 



As shown in the phase plane diagram, the system of concern does 
not converge but circles the origin. Since the equation was that of 
an undamped oscillator, its response should be as shown or simple har- 
monic motion. 



3-2 Application of Phase Plane Analysis to the Problem of Transient 
Response of the Delay-Lock Discriminator. 

The differential equation describing the transient response of the 
Digital Delay-Lock System will now be derived. Figure 3-2 shows a simpli- 
fied version of the system under consideration. 



Sn,(i+f*Y) 

-5« 




Fig. 3-2 Simplified block diagram of the Digital 
Delay- Lock Discriminator. 

The loop filter for this system is a simple low pass RC network as 
shown in Fig. 3-3 along with the filter transfer ration. 




Fig. 3-3 
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The system equation can be written in terms of the estimate of the Delay 



derivat ive 
» 

A 



t . 



and the Discriminator Characteristic 






such that 



* where is normalized to be ^ ^ 



pee) 



Expressed in terms of operator notation, the equation becomes: 



pf ' Rp)- [p(£)] 



where p denotes a time derivative operator J 



Expanding the above equation gives: 

P]t= D(e) + p[^D( 0 ] 



\ i'i:)- Vo ' i To 
6(t) =TCt)-f ("t) 



then 



or 



T = D(t) tl/F D'(0-e 
'£ ’ P(0+l/?D'fe)'€ 
let ^ slope 

4 

-T6 PCfc) + VT D(t)-6 
DC6) 



£ = Vo- T 

T. 



£ =-t 



' c = 



r+i/Fd'a) 



For a finite loop gain system, the filter transfer function was 



choosen as 




)-fVT p 



and in terms of operator notation, the 



system equation becomes: 

p T tr T( p) • DCt) 

[pV.iF]t= 

f+.iT = D(fc)^)/ 7 D‘(y-E 

+A(v-'£)-ffD'(0-e * D(e) 

or ; _ D (0 - 0-1 /. 

^ ~ V+o-i ti/? D“(0 

These then are the system equations which when used to construct the 



phase plane, will describe the transient response of the system. In 
order to construct a phase plane that will be of use, it is necessary 



to evaluate the equation for numerous combinations of slope and 

delay error ^ . The next topic to be considered describes the use of 
a high speed computer to construct the required phase planes. 



3-3 Construction of Phase-Planes by Digital Methods. 

As previously mentioned, construction of the phase plane for a 
given system can be a time consuming task of point by point evaluation 
to plot the desired isoclines. When the system is a non-linear one, this 
task becomes even more difficult as there may be no obvious symetry to the 
phase plane. Therefore, the problem of evaluating the required isoclines 
to produce transient response trajectories, is simplified when accomplish- 
ed on a high speed computer. Before describing the results of the com- 
puter program, it is first necessary to further define the Discriminator 
Characteristics of interest. 

As mentioned in section 2, the digital system has a saw tooth character- 
istic as shown in Fig. 3-4, and is mathematically defined as follows: 



D(0 = 



— 6 •^'oy 

^ for 

- 6+2 

o else 



- 2^6 ^-1 
-U<€ < I 
K 6 
U>V\erd 







Fig. 3-4 Saw tooth Discriminator Characteristic 
The Gaussian Discriminator Characteristic as previously shown in Fig. 

2-3 is also of interest in evaluating the basic delay-lock system. This 
function is mathematically described as follows: 
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The computer programs required to evaluate the isoclines are de- 
tailed in Appendix A; it is only necessary tc say here that many 
hundreds of points were required to complete the phase planes. The 
plot of the isoclines for the saw tooth characteristic is shown in 
Fig. 3-5. It is apparent that the isoclines of magnitude greater than 
5 are concentrated near the ^ axis and so could not be plotted. A 
sample trajectory is sketched on the phase plane and shows the lock-on re- 
sponse for one set of initial conditions. It is apparent that unless a 
greater number of isoclines is plotted, the accuracy of the trajectory is 
to be deeply suspect. This however is the normal method of solution by 
phase plane techniques. Clearly, this method leaves much to be desired 
and the next section will show two levels of improvement by evaluating 
not only the phase plane but also any specified trajectory, with the 
digital computer. Fig. 3-6 shows the phase plane constructed for the 
Gaussian Characteristic. 
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Transient Analysis By Sequential Programming , 



As previously mentioned, the transient response of a given system 
may be obtained by constructing a phase plane associated with the 
system and then sketching the required trajectories by faring the 

curve through the known isoclines. This technique has been recently 

• /3/ 

applied to a system known as the Phase-Lock Discriminator ; however, 

this method is at best adequate in that the errors involved in sketch- 
ing are acqumulative, and the final value for lock-on lock critical 
trajectories may be meaningless. Certain problems can be adapted to 
solution by analog computer methods^^^, but here again one has to worry 
about amplifier drift during solution producing flucuating errors. For 
the solution of the non-linear equation of the delay-lock, analog methods 
are not practicable, and some other method is desirable. 

This section will present a method of digital analysis by computing 
the desired trajectories with the aid of a high speed general purpose 
digital computer, the CDC 1604. 

Two levels of improvement will be covered; the fixed increment 
approach, and the method of Dynamic or Sequential Programming which 
limits the error to any fixed amount. Actual computer results for the 
transient response of the Delay-Lock Discriminator will be presented. 
Although time is not directly available in the phase plane diagram, it 
is inherently involved and can be obtained by graphical procedures^^^ . 
However, it is available during the computer solution and this method will 
also be covered. 

The Flow Diagram and Fortran Source Language Programs are contained 
in Appendix B. 

4-1 Solution of the Differential Equation by Fixed Increment Method. 

The equation to be solved by phase plane methods as described in 



29 



• DCO 

section 3 is given by: ^ 



Rearranging to solve for the slope gives: 

This is the main value to be computed as the solution of the trajectory 

developes* Figure 4-1 shows the basic method of approach. Starting with 

a known point in the phase plane representing the initial conditions, 

. 1^0 ) 

the slope. Qq at this point is evaluated. Using a fixed increment of ^ , 

a£ , the next value of velocity error ^ is computed in a linear 

• • * 

fashion such that and the value of delay at this 

point is determined by . This process is continued 

until the trajectory reaches the stable point (in this case the origin) 
where the system is then in lock and tracking, or until the trajectory 
reaches a value of delay error ^ equal to 2 and the system fails to lock 



in. 




Fig. 4-1 Computer Approach to Evaluation of Trajectory. 

In the method just described, the increment value was chosen to be 
= 0.01. Although this method produces satisfactory sets of tra- 
jectories, there is a major difficulty involved, in that when the tra- 
jectory approaches the ^ axis, the slopes in the phase plane are 
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approaching infinity and even with a very small inclement^ the new value 

of ^ will be made to go too far negative and not be the true trajectory. 

An example of this is shown below and taken from a sample computed trajectory. 

When the value of first went negative, the routine would set the new 

• • • 
value of £: , , equal to the negative of the old value and 

then resume the program. 



Value of ^ 



Value of € Slope 



1.11000 0.29645 -9.07462 

1.12000 0.20569 -13.75419 

1.13000 0.06815 

switch over occurs as next value 
of ^ is negative 

1.13000 -0.06815 122.299 

1.12000 -1.29114 3.83068 

Note this large increase in ^ in 
just one value of increment^ At# 



Figure 4-2 shows the results of this first approach to the evalua- 
tion of the transient response of the system. Several values of initial 
conditions were plotted to describe the lock-on region of interest. 



4-2 Sequential Programming Applied to Phase Plane Analysis. 

Sequential or Dynamic Programming is a basic adaptation of the pro- 
cess developed by Bellman^^^^. It is generally described by consider- 
ing the computer as a adaptive device capable of changing system para- 
meters internally as the solution of a given problem develops. In the 
case of the computer solution of the phase plane trajectories, this means 
not operating with a fixed increment, but rather adapting a variable in- 
crement to follow the solution as it develops. 

Instead of the fixed increment as previously covered, consider a 
variable increment a€k , which will take on one value as the slope 
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in the phase plane becomes zero, and decreases to zero as the slope 
approaches infinity. It was decided that the variable increment in 



this problem would be defined as follows: 




This allows an increment to take on values from 0.02 to zero. When the 
trajectory passes through slope areas approaching infinity, the trajectory 
starts to move almost vertically until this high slope region is passed. 

A comparison of this method with that of the fixed increment type is 
shown in Fig. 4-3. The method of sequential programming then allows the 
trajectory to recover from these high slope regions without being overly 
influenced by them. 



Figure 4-4 shows the results of variable inclement method for deter- 
mining transient response of the Delay-Lock System. These curves are based 
on a saw-tooth discriminator characteristic and the system having infinite 
loop gain. Various initial velocity error conditions were plotted to 
demonstrate the lock-on characteristics of the system. Considering the 
trajectory describing an initial velocity error of unity, the curve starts 
to increase in velocity error while the system is reducing the delay error 
and correcting the estimate of signal delay. The trajectory starts its 
spiral towards the singular point which corresponds to a lock-on and track 




N 






\ 



Fig. 4-3 Comparison of fixed and variable 
increment trajectories. 
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condition. In this case the system does lock on and so will tolerate 
a velocity error of this magnitude. However, in the case of an initial 
velocity error of 2.3, the system attempts to lock-in and track but 
does not quite make it and assumes a no- lock condition. The system 
will then remain out of lock until the estimate of the received signal 
is again in the discriminator region. Fig. 4-5 shows the transient 
response for the same system discriminator characteristic but concerns 
the finite loop gain case. The lock-no lock critical trajectory occurs 
for the same initial velocity error but in this latter case, the system 
comes to rest with a finite delay error proportional to the initial 
velocity error. Since a trajectory may be started any where in the 
phase plane, the computer curves offer a second analysis of the operating 
system. This concerns the case where the system is in lock and tracking 
a target with a constant velocity. If the target accelerates (or de- 
celerates), then the target assumes a new velocity and the system is 
temporarily out of track and must lock on in the manner described above. 

The transient response for the delay- lock system with a Gaussian 
Discriminator Characteristic and infinite and finite loop gains is shown 
in Fig. 4-6 and 4-7. The phase plane curves are explained in the manner 
previously mentioned. Although the computer data was hand drawn, the 
curves could easily be constructed by a high speed digital plotter to 
give as many trajectories as desired. 

The error in the trajectory as it passes through the region of large 
valued slopes may be approximated in the following manner; 



Assuming a parabolic trajectory when crossing through the infinite slope 
area allows the following simplifications: 



The next value of t mav be expanded in a Taylor's series to be: 
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4-3 Computation of the Lock-On Time Response. 

As previously mentioned, time response is not directly available 
for the phase plane diagram, but may be determined by graphical tech- 
niques. However, this procedure is again a tiresome method which is not 
very accurate. It is possible to compute the time response for a given 
trajectory while the computer routine is in the process of constructing 
the trajectory. Fig. 4-9 shows the typical delay error versus time response 
for several initial velocity conditions along with a typical velocity error 
versus time response. As is apparent, the delay error term starts 

at a value of -2 and oscillates in a damped fashion converging to zero 
value delay error when the system has attained Lock-On. These transient 
responses were observed on a analog plotter for the Digital System and 
were of the same general appearance as those predicted. Fig. 4-8 shows 
an observed transient response for the approximate velocity error function 
when the system was in the process of lock-on. In order to show the method 
of computing transient times it is necessary to expand the delay error 
function in a Taylor's Series. Therefore, a)= + 

zi 

and the increment of time required for the trajectory to pass from 
to ^ is given by: — | ^ 



1 
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It is only necessary to sum the time increments to find the total 



time expired at any given point on the trajectory. The small error 

2 * * 

term At € was neglected. 




Fig. 4-8 Observed wave-form of velocity error 
response of Digital Delay-Lock System 
during lock-on. 
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5 . Conclusions, 



The Digital Delay-Lock Discriminator with its pseudo- random signal 
characteristics offers a method of tracking targets at great ranges and 
high velocities, while operating with moderate power requirements. The 
experimental system is now operating at the Lockheed Missle and Space 
Company and is under consideration by NASA as a proposed tracking and 
rondevous system for maned and unmaned satelites. The system has other 
possibilities such as an underwater tracking system. Since the system 
has the capability to track targets when the Signal to Noise ratio is 
in the order of -30db, the power of the transmitted signal could be 
maintained at a level below that of the ambient noise surrounding the 
target, and as such, would not give the target an indication it was 
under survailence. 

In order that the system be optimized, it was necessary to deter- 
mine the transient response of a normalized system during the lock-on 
process. Phase-plane analysis methods were then utilized to obtain 
this response. With out the aid of a high speed digital computer, this 
analysis would not have been easily accomplished. The data obtained 
from the sequential programming methods demonstrate that this phase 
plane approach to the analysis problem is applicable to a large number 
of non-linear systems. 
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APPENDIX A 



This Appendix contains the Flow Diagram and Fortran Source 
Programs required to evaluate the isoclines (lines of constant slope) 
of the phase planes which represent the solution to the differential 



equation: 



-e = 



])ce) 






In this problem, it is necessary to fix a value of slope , 

and evaluate the above equation for £ with €. taking on the range 

; an increment of 0.01 was choosen for 
Values of slope ^ used were: + ^ 0.2, 0.4, 0.6, 0.8, 1.0, 2.0, 3.0, 

4.0, 5.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0 
It was thought that this many values of slope would be required to 
properly bound the phase plane. As it turned out, only slopes in the 
region 10 were useful since those of higher value could not be 

plotted on a useable phase plane. The results of these programs are 
shown in Section 3, Figures 3~5 and 3~6. 
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PROGRAM DELAY 1 

THIS PROGRAM PRODUCES REQUIRED 150 CLINES 
FOR SOLUTION OF THE DELAY-LOCK DISCRIMINATOR 

PRIN'l 22 
V=0.04 
W=K-2 
V=W<V 

PRINT 33.V I 

DO 14 1=1,36 

N=5 

M=9 

L= 18 

NM=23 

ML =27 

LN=36 

IF(SLOPE-50.0)51,61,51 
61 SLOPE=0.0 
51 IF(I-N)3,3,1 

1 IF(I-M)4,4,2 

2 IF( 1-L)5,5,23 

3 SL0PE=SL0PE+0.2 
GO TO 6 

4 SL0PE=SL0PE+1.0 
GO TO 6 

5 SL0PE=SL0PE+5.0 
GO TO 6 

23 IF( I-NM)26,26,24 

24 IF( I-ML)27,27,25 

25 IFI I-LN}28,28, 15 

26 SL0PE=SL0PE-0.2 
GO TO 6 

27 SLOPE=SLOPE-1.0 
GO TO 6 

28 SLOPE=SL0PE-5.0 

6 PRINT 55, SLOPE 
DO 13 J=l,81 
H=J 

- £=-4.1+0.1*H 
EI=E*E 

C=2. 7182818284 
P01=E»( l.O/C**El) 

PD2=( 1.0-2.0*E1 )*{ 1.0/C**E1) 

12 Y=l PD1+SQRTF(2.0)*PD2*V) /I SLOPE+SQRTF (2.0)*PD2 ) 

13 PRINT 44,E,Y 

14 CONTINUE 
PAUSE 

50 CONTINUE 

15 STOP 

22 FORMAT(9H VELOCITY, 3X, 6H SLOPE, 3X,8H E-VALUE, 15X,8H Y-VALUE,//) 
33 F0RMAT(F8.4, /) 

44 FORMAT(25X,F4.2, 16X,F9.5,/) 

55 FORMAT! 12X,F8. 4,/) 

END 

END 
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PROGRAM DELAY 2 

C THIS PROGRAM PRODUCES RECUIRHD 150 CLINES 

C FOR SOLUTION OF THE DELAY-LOCK DISCRIMINATOR 

DO .50 K=2, 12,2 



um 



61 

51 

1 

2 

•a 



23 

24, 

25' 

26 

27 

28 
6 



8 

9 

10 

1 1 
12 
.13 

14 

50 

15 
22 
33 
44 
55 



V=0.04 
W=K-2 
V=W*V 

PRINT 33, V 
DO 14 I=U36 
N=5 
M = 9 
L= 18 
NM=23 
ML = 27 
LN=36 

IF{SL0PE-50.0)51,61,51 

SLOPE=0.0 

IF(I-N)3,3,1 ' 

IF( I-M)4,4,2 
IF( I-L)5,5,23 
SL0PE=SL0PE+0.2 
GO TO 6 

SLOPE=SLOPE+ 1 .0 
GO TO 6 

SL0PE=SL0PE+5.0 
GO TO 6 

IF( I-NM)26,26,24 
IT=( I-ML)27,27,25 
I I-LN)28,28, 15 
SL0P£=SLOPE-0.2 ' 

GO TO 6 

SL0P6=SL0PE-1.0 
GO TO 6 

SLOPE-SLOPE-5. 0 
PRINT 55, SLOPE 
DO 13 J=U 
DO 13 J=l,38 
H=J 

E=-2,0 +0. 1*H 
TEMP=ABSF( E ) 

IF(T6MP-1.0)7,7,8 

P01=E 

GO TO 9 

PDl = Ei- 2.0*E/TEMP 
IF ITEMP-1.0) 10, 10, 11 
PD2- 1.0 
GO TO 12 
PD2 1.0 

Y=( PD 1+SQRTFI 2.0)*P02*V) / ( SL0PE+ SQRTF ( 2 . 0 ) * PD2 ) 

PRINT 44,E,Y 

CONTINUE 

PAUSE 

CONTINUE 

STOP 

F0RMATI9H VELOC I TY , 3X , 6H SLOPE, 3X,8H E- VALUE , 1 5X , 8H Y-VALUE,//) 
FORMAT (F8. 4,/) 

FORMAT (25X,F4. 2, 16X,F9.5,/) 

FORMAT! 12X,F8. 4,/) 

END 

END ' ‘ 
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APPENDIX B 



This appendix contains the Flow Diagrams and Fortran Source 
Programs required to generate the desired trajectories describing the 
lock-on behavior of the system. The results of these programs are 
shown in Section 4, Figures 4-2 through 4-7. 
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1 

READ ],X2 

2 READ 1, Y 
X=X2 

PRINT 40 

3 E1=X*X 
0=2.7182818284 
PD1 = X*( 1.0/C**E1 ) 

PD2=t 1.0-2.0*E1 )*( 1.0/C**E1) 

A=SCRTF12.0) 

BET/s = (PDl + (PD2*Y*A) )/(-Y) 

TEMP=ABSF(BETA) 

IFITEMP-0.5)31,31,32 

31 DEL1A=0.1 
GO TO 33 

32 0ELTA=0.01 

33 Y1=Y+BETA«DELTA 
X1=X+DELTA 
IF(Xl-2.0)34,34, 13 

34 IF(Y1)5.5,4 

4 PRINT 50, XI, Yl, BETA 
Y=Y \ 

X=X 1 
GO TO 3 

5 IF(DELTA-0.01 )7,7,6 

6 DELTA=0.01 
GO TO 3 

7 D6LTA=-0.01 

IF( X-0. 1 ) 12, 12,8 

8 Y=-0.025 
81 E1=X*X 

0=2.718281284 
P01=X*( 1.0/0**E1 ) 

P02=( 1.0-2.0*E1 )*( 1.0/0»»l:l ) 
BETA=(PDl+(P02*Y*A))/(-Y ) 

Y1=Y+BETA*DELTA 
X1=X+DELTA 
IF(Y1)9, 10, 10 

9 PRINT 50,X 1,Y1 ,BETA 
Y=Y I 

X=X 1 

GO TO 81 

10 IF(X 1+0. 1 ) 1 1, 12, 12 

1 1 Y=0.025 
'^ELTA=0.01 

GO TO 3 
,12 PRINT 70 
' PAUSE 
* GO TO 2 
13 PRINT 60 
PAUSE 
GO TO 2 

40 F0RMAT(8H X-VALUE , SX, 8H Y~VALUE,5X,5H BETA,//) 
50 F0RNAT(F8.5,5X,F8.5,5X,F8.5,/) 

60 FORMAT! IIH NO LOOK-ON,//) 

70 FORMAT (8H LOOK-ON,//) 

END 

END 



^2 




^ 0 , 5 ) 



X FORRM ^ 

2 READ J,Y' 

X=X2 

PRINT 40 

3 TEMP=ABSF{X) 

IF(TEMP-].0)4,4,5 

4 PD1=X 
P02=1.0 
GO TO 6 

5 PD1=X+2.0*X/TEMP 
PD2=-1.0 

6 A=SQRTF(2.0) 

BETA={PDl+(PD2*Y*A))/(-Y» 

TEMPO=ABSF( BETA) 

DELTA=0.02/( l.O+TEMPO) 

Y)=Y+3ETA*DELTA 
X1=X+DELTA 
IF(Xl-2.0)7,7, 19 

7 IF(Y1)9.9,8 

8 PRINT 50, XI, Y1 .BETA 
Y=Y I 
X=X1 
GO TO 3 

9 CONTINUE 
IFIX-0. I ) 18, 18, 10 

10 Y=-Y 

11 TEMP=ABSF(X) 

IF(TEMP-l.O) 12, 12, 13 

12 P01=X 
PD2=1.0 
GO TO 14 

13 PI)l = Xi-2.0*X/TEMP 
P02=-1.0 

14 BETA-^ PD1+(PD2*Y*A) )/(-Y) 

TEMPO"ABSF(BETA) 

- DELTA=—0.02/( l.O+TEMPO) 

Y1=Y+BETA*DELTA 
X1=X+0ELTA 
IF(Yl) 15, 16, 16 

15 PRINT 50, XI, Yl, BETA 
Y=Y1 
X=X1 

GO TO 11 

16 IFIXI^-O. 1) 17, 18, 18 

17 y=-Y , ’ 

GO TO 3 

18 PRINT 70 

PAUSE •, ^ 

GO TO 2 ! 

19 PRINT 60 
PAUSE 2 
GO TO 2 

40 FORMATION X-VALUE, 5X, 8H Y-VALUE, 5X,5H BETA,//) 
50 F0RMAT(F8.5,5X,F8.5,5X,F8.5,/) 

60 FORMAT! 1 IN NO-LOCK-ON,//) , 

70 F0RMATI8H LOCK-ON,//) f 

END ! ! . 

END i ; i 
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PROGRAM TJ^AKON 
1 FORMAK F10.5) 



READ 



vX2 

«.y > 



I 



3 



4 

5 



X=X2 

PRINT 40 
FLAG=0..0 
E1=X*X 

C=2.71«2818284 
PD1 = X<»1 1.0/C**E1) 

P02={ 3..0-2.0*El )*( 1.0/C**6 1 ) 

A=SQRTF(2.0) 

BETA=( ( PDl+( 0. 1*V) )♦{¥*( A*P02+0. !)))/(-¥) 
TEMPO=ABSF(BETA) 

DELTA=0.02/( l.O+TEMPO) 

Y1=Y+BETA*DELTA 

X1=X+DELTA 1 

IF(Xl-2.0)4,4, 12 

IF(Y1)6.6,5 

PRINT 5.6, XI, Y1, BETA, V 

Y=YI 

X=X 1 

GO TO 3 



6 IF(FLAG.)7,7, 1 1 

7 FLAG=1.0 
Y=-Y 

8 E1=X*X 

PD1=X*( U0/C**E1) 

• PD2=( U0-2.0*E1)*( 1.0/C*»E1 ) 

BETA={(PDl+(0. 1*V) )+(Y*{ A*PD2+0. 1) )) /(-Y) 

TEMPO=ABSF(BETA) 

DELTA=-0.02/( l.O+TEMPO) 

Y1=Y+B6TA*DELTA 
X1=X+DELTA 
IFIY1 )9. 10, 10 

9 PRINT 56,X1,Y1,BETA,V 
Y=Y1 

X=X1 

GQ^O 8 ' 

10 Y^ 

GO TO 3 

1 1 PRINT 70 
; PAUSE 

GO TO 2 

12 PRINT 60 
PAUSE 2 
GO TO 2 

40 F0RMAT(8H X-VALUE , 5X, 8H Y-VALUE ,5X,5H B£TA,5X,8H V-VALUE,//) 
50 FORMAT! F8.5,5X,F8.5,5X,F8-5,7X,F 10-5,/) 

60 FORMAT! IIH NO-LOCK-ON,//) 

70 F0RMAT!8H LOCK-ON,//) 

END 

END 
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9 

10 

1 1 

12 

13 

14 



15 

16 
18 

19 



40 

SO 

60 

70 



PROGRAM TRAKIN 
FORMATC F 10.5) 

Rise 1:?*^ I 

V=-Y 1 

X=X2 i 

PRINT 40 
FLAG=0.0 

TEMP=ABSF(X) < 

IF(TEMP-1.0)4,4, 5 
PD 1 = X 
PD2=1.0 
,.G0 TO 6 

PD1=X+2.0*X/TEMP 
PD2=-1»0 
A=SQRTF(2.0) 

BETA=( (PDl + (0. 1*V) ) + (Y*{ A*P02+0. 1 m /(-Y) 

TEMPO=ABSF(BETA) 

OELTA=0.02/( l.O+TEMPO) 

Y1«Y+8ETA»DELTA 
Xli^x + DELTA 
IF(Xl-2.0)7,7, 19 
IF( Y1 )9.9,8 
PRINT 50,X1,Y1,BETA,V 
Y=Y I 

X=X 1 i 

GO TO 3 

IF( FLAG) 10, 10, 18 
FLAG=KO 
Y=-Y 

TEMP=ABSF(X) 

1F( TEMP-1.0) 12, 12, 13 
PD1 = X 
PD2=1.0 
GO TO 14 

PD1=X+2.0»X/TEMP f 

PD2=-1.0 

BETA=( ( PDl+( 0. 1*V) )+(Y*( A*PD2+0. 1) )) / (-Y) 
tEMPO=ABSF( BETA) 

DELTA=-0.02/ ( 1 .O+TEMPO) 

Y1=Y+BETA*DELTA 
X1=X+DELTA ! 

IF( Y1 ) 15, 16, 16 i 

PRINT 50, XI , Y 1 ,BETA, V 

Y=Y 1 ; . 

X=X1 

GO TO 1 1 ! 

Y=-Y i 

GO TO 3 ; 

PRINT 70 I 

PAUSE 

GO TO 2 

PRINT 60 

PAUSE 2 

GO TO 2 

F0RMATJ8H X-VALUE, 5X, 8H Y- V A LUE , 5X ,5 H BETA,5X,8H V-VALUE,//) 
FORMATC F8.5,5X,F8.5,5X,FE.5,7X,F 10.5 ,/) 

FORMATC 1 IH NO-LOCK-ON,// ) , 

FORMAT C8H LOCK-ON,//) j 

END i 

END i 
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